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1 INTRODUCTION 



ABSTRACT 

We present the results of a deep K-band imaging study which reveals the host galaxies 
around a sample of luminous radio-quiet quasars. The -ftT-band images, obtained at 
UKIRT, are of sufficient quality to allow accurate modelling of the underlying host 
galaxy. Initially, the basic structure of the hosts is revealed using a modified Clean 
deconvolution routine optimised for this analysis. 2 of the 14 quasars are shown to 
have host galaxies with violently disturbed morphologies which cannot be modelled 
by smooth elliptical profiles. For the remainder of our sample, 2D models of the host 
and nuclear component are fitted to the images using the \ 2 statistic to determine 
goodness of fit. Host galaxies are detected around all of the quasars. The reliability 
of the modelling is extensively tested, and we find the host luminosity to be well con- 
strained for 9 quasars. The derived average if-band absolute iC-corrected host galaxy 
magnitude for these luminous radio-quiet quasars is (Mk) = —25.15 ± 0.04, slightly 
more luminous than an L* galaxy. The spread of derived host galaxy luminosities is 
small, although the spread of nuclear-to-host ratios is not. These host luminosities 
are shown to be comparable to those derived from samples of quasars of lower total 
luminosity and we conclude that there is no correlation between host and nuclear 
luminosity for these quasars. Nuclear-to-host ratios break the lower limit previously 
suggested from studies of lower nuclear luminosity quasars and Seyfert galaxies. Mor- 
phologies are less certain but, on the scales probed by these images, some hosts appear 
to be dominated by spheroids but others appear to have disk-dominated profiles. 
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suggested to be super -massive black holes potential l y capa- 



Models of the cosmological evolution of quasars often use 
galaxy mergers as the primary mechanism for quasar acti- 
vation and require the mass of the structure within which a 



ble of powering A GN ( Kormendy fc Richstone 1995 ; Magor 



rian et al. 1998). These studies have also found evidence 



quas a r is formed as a basic p arameter ( Efstathiou fc Ree, 1 



198S ; Haehnelt & Rees 1995 



Percival & Miller 



1999p . 



One 

step towards testing hypotheses about quasar initiation is 
to answer the question: Is a quasar's luminosity correlated 
with the luminosity of the structure within which it formed? 



Such a correlation has been shown to exist for low red- 
shift (0 < z < U.3) Seyterts and quasars with luminosities 
Mv ^ —25 in that there appears to be a lower limit to 
t he host luminosity wh i ch increases with quasar luminosity 
(McLeod fc Rieke 1995; McLeod, Rieke & Storrie-Lombardi 



for a correlation between the mass of the compact object 
and the luminosity of the spheroidal component of the host. 
Assuming a link between nuclear luminosity and black hole 
mass, the average nuclear luminosity emitted by low redshift 
quasars is expected to increase with host spheroidal lumi- 
nosity. In light of this prediction there has been a re surgence 
of interest in host galaxy studies and recent work ( McLure 



et al. 199E ) has found weak evidence for a correl ation in ac- 

However 



199£). However this limit is poorly defined, particularly for 
high luminosity quasars when the strong nuclear component 
makes it increasingly more difficult to find low luminosity 
hosts. 

Recent work has shown that the majority of nearby 
galaxies have massive dark objects in their cores, which are 



cord with the relations of Magorrian et al. (1995 
this correlation relies on spheroid/disk decomposition for 
two quasars with low nuclear luminosities, and only a small 
number of luminous radio-quiet quasars were observed. 

Host galaxy properties of AGN are known to be corre- 
lated with the radio power: radio galaxies tend to be large 
spheroidal galaxies, while disk galaxies tend to be radio- 
quiet. Recent evidence suggests that the hosts of radio-loud 
quasars are also predominantly massiv e spheroidal galax- 
ies regardless of the nuclear luminosity ( Dunlop et al. 199E ; 
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Taylor et al. 1996; McLure et al. 1999). However, studies 
of radio-quiet quasars with luminosities My > —25 have 
shown that the hosts can be dominated by either disk-like 
or spheroidal components or can be complex systems of gray 



itation ally interacting components ( Taylor et al. 1996 



call et c I. 1997; Boyce et al. 1995 



McLure et al. 



Bah- 



There 

is therefore strong justification for studies to see if these lu- 
minosity and morphological trends extend to the hosts of 
more luminous (Mv ^ —25) radio-quiet quasars. 

There have been many recent detections of host galaxies 



in the optica l thanks to results f rom HST (Hooper, Impcy 



& Foltz 



1997 



Bahcall et al. 1997; Boyce et al. 1998 



McLure 



et al. 1 $9|)7 which add to our knowledge from ground-based 



studies Jpalkan et al. 1984a|; |Malkan 1984b| ; [Hutchings & 
Neff 19p2|; |Veron-Cetty fc Woltjcr 199C^ 7Tiowever, quasar 
hosts often appear significantly disturbed, as if by inter- 
action or merger which can lead to strong bursts of star- 
formation and significant extended line and blue continuum 
emission at optical wavelengths which are not indicative of 
the mass of the underlying host. The nuclear-to-host light 
ratio in the optical is also typically higher than at longer 
wavebands. 

These problems can be circumvented by observing in 
the infrared where the contrast between host and nuclear 
component is improved and the emission associated with 
star bursting activity is largely absent: the K magnitude is 
a better m easurement of the long-li ved stellar populations 
in the host (Bruzual & Chariot 1993). Previous observations 
in the infra-red have been successfully used to determine 
quasar host galaxy luminosities and morphologies (McLeod 
& Rieke 1994a;b; Dunlop et al. 1993; Taylor et al. 1996). 
However, recent advances in telescope design, in particu- 
lar the advent of adaptive optics systems such as the tip-tilt 
system on UKIRT produce clearer images of quasars and en- 
able accurate point spread functions (psfs) to be more read- 
ily obtained as differences between successive observations 
are reduced. Such advances coupled with improved analysis 
techniques mean we are now able to reveal the host galaxies 
of luminous quasars with Mv < —25 using infra-red obser- 
vations. 

In order to obtain enough luminous radio-quiet quasars 
our sample was forced to cover redshifts 0.26 < z < 0.46. 
At such redshifts, with typical seeing, the structure of the 
host galaxy is hidden in the wings of the psf from the nu- 
clear component. There are two main ways of proceeding: 
either the psf can be deconvolved from the quasar light to 
directly observe the host galaxy, or known galaxy profiles 
can be used to model the hosts, a nuclear component can be 
added in and the profiles can be fitted to the data. Because 
the host galaxies sometimes have disturbed morphologies in- 
dicative of violent mergers, it is difficult to assume a form 
for the galaxy. However without such modelling, it is not 
easy to determine the contribution of the host to the light 
from the centre of the quasar and deconvolution routines 
tend to produce biased solutions which may alter important 
features. 

For the analysis of our quasar sample, an approach 
is adopted which uses both methods. Initially (Section [s]) 
the images were restored using a deconvolution algorithm, 
based on the Clean algorithm (Hogbom 1974), devel oped 



for thi s problem, w hich will be described elsewhere (Per 



extent to which the 'nebulosity' around the point source is 
disturbed. Deconvolution of the light from two of our quasars 
reveals violently disturbed host galaxies indicative of close 
merger events. In the remainder of our sample, the non- 
nuclear light is more uniformly distributed around the centre 
of the quasar. We should note that the resolution provided 
by this deconvolution technique is probably not sufficient to 
reveal evidence for weak mergers, where the host galaxy is 
only slightly disturbed. 

Where the image-restoration routine revealed approxi- 
mate elliptical symmetry in the non-nuclear component, 2D 
galaxy profiles were fitted to the hosts. Analysis of non- 
interacting, low redshift galaxies has shown that an empiri- 
cal fit to both disk and spheroidal systems is given by: 



/i = /i D exp 



i/i3 



(1) 



where /i is the average surface brightness in concentric ellip- 
tical annuli around the core, and r is the geometric average 
of the semi-major and semi-minor axes. 

Model images were carefully created using this profile 
and were tested against the data using the \ 2 statistic to de- 
termine goodness of fit. Five host parameters were required, 
the half-light radius, integrated luminosity, axial ratio, angle 
on the sky, and the power-law parameter of the galaxy f3, 
as well as the nuclear-to-host ratio. Section [| describes the 
modelling procedure in detail, and in Section ^ the best-fit 
parameters are presented for the host galaxies. 

Much previous work has produced ambiguous results 
because of a lack of error analysis and insufficient testing of 
the modelling. A detailed analysis of the reliability of the 
2D modelling method used in this paper has therefore been 
undertaken and is presented in Section ^. Although hosts 
are detected in all of our sample, the upper limit of the 
host luminosity is only usefully constrained for 9 of the 12 
quasars modelled (see Section 5.1). Similar analysis of the 
best-fit (3 parameter which determines the morphology of 
the host reveals that this parameter is, unsurprisingly, more 
poorly constrained than the luminosity. However, we have 
created Monte-Carlo simulations of images with the same 
signal-to-noise as the original images (Sections | By 
analysing these images using exactly the same procedure as 
for the original data we find that it is possible to distinguish 
between disk and spheroidal structure. 

Unless stated otherwise we have adopted a flat, A = 
cosmological model with Ho = 50kms _1 Mpc _1 and have 
converted previously published data to this cosmology for 
ease of comparison. 



2 THE SAMPLE AND OBSERVATIONS 

We have selected 13 luminous (Mv < —25.0) quasars and 
one less luminous quasar within the redshift range 0.26 < 
z < 0.46. The quasars were checked for rad io loudness us- 
ing the NVSS survey (|Condon et al. 199Sj). Three of the 



cival fc Miller 200C ). This routine was used to reveal the 



14 quasars were detected at 1.4 GHz in this survey (see Ta- 
ble |l|) but their flux densities are all < 10 24 W Hz^Sr" 1 
and they are considered part of the radio-quiet population. 
Three quasars, 0137-010, 0316-346 and 2233+134 were ob- 
served at UKIRT before the tip-tilt system was operational 
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quasar J2000 coords V z My 1.4 GHz Flux density Observing run Integration time 

/WHz-'Sr" 1 /s 



0043+039 


00 45 47.3 


+04 10 22.5 


16.0 


0.384 


-26.0 


- 


09/1997 


2800 


0137-010 


01 40 17.0 


-00 50 03.0 


16.4 


0.335 


-25.3 


1.46X10 23 


09/1996 


11300 


0244-012 


02 46 51.8 


-00 59 32.3 


16.5 


0.467 


-25.9 


- 


09/1997 


10300 


0316-346 


03 18 06.5 


-34 26 37.1 


15.1 


0.265 


-26.0 


- 


09/1996 


6400 


0956-073 


09 59 16.7 


-07 35 18.9 


16.5 


0.327 


—25.1 




05/1998 


8000 


1214+180 


12 16 49.1 


+17 48 04.1 


16.7 


0.374 


-25.2 




05/1998 


7590 


1216+069 


12 19 20.9 


+06 38 38.4 


15.7 


0.334 


-26.0 




05/1998 


7200 


1354+213 


13 56 32.9 


+21 03 51.2 


15.9 


0.300 


-25.5 




05/1998 


8600 


1543+489 


15 45 30.2 


+48 46 08.9 


16.1 


0.400 


-26.0 




05/1998 


14300 


1636+384 


16 38 17.6 


+38 22 49.0 


17.0 


0.360 


-24.8 




05/1998 


7100 


1700+518 


17 01 24.9 


+51 49 20.4 


15.1 


0.290 


-26.2 


7.18X10 23 


09/1997 


3600 


2112+059 


21 14 52.6 


+06 07 42.5 


15.5 


0.466 


-26.7 


2.95X10 23 


09/1997 


15100 


2233+134 


22 36 07.7 


+13 43 55.0 


16.7 


0.325 


-26.9 




09/1996 


6400 


2245+004 


22 47 41.6 


+00 54 57.3 


18.5 


0.364 


-23.4 




09/1997 


11100 



Table 1. The sample of 14 quasars observed. My was calculated for each quasar from apparent magnitudes given in the catalogue 
of Hewitt & Burbidgc ( |l993j ) assuming no ftf-correction. Radio fluxes were determined using the NVSS 1.4 GHz survey. Only 3 of the 
quasars were detected in this survey and their radio fluxes are below the radio-quiet/loud cutoff. We observed this sample in 3 observing 
runs at UKIRT. For the 09/1996 run, the tip-tilt system was not operational and the psf stars are not expected to be as well matched 
to the quasars as for the other runs (see Section ^) . Note that coordinates given were determined using the original finding charts and 
the Digitised Sky Survey and may be different from those in the catalogue of Hewitt &; Burbidge (1993) which are often inaccurate. 



and so these data are not of the same quality as those from 
subsequent runs. 

Of the 13 luminous quasars selected, three, 0956—073, 
1214+186 and 1636+384 have not had previous attempts to 
measure host magnitudes and morphologies. It is difficult 
to assess the significance of claimed host detections for the 
other quasars and the associated parameters calculated be- 
cause of the lack of error analysis which abounds in this field, 
and the great potential for systematic errors caused by the 
requirement for accurate psf measurements. However, indi- 
vidual results from thes e stu dies are compared to the results 
of this paper in Section 



8.4 



The observations were all taken using the 256 x 256 pixel 
InSb array camera IRCAM 3 on the 3.9 m UK Infrared Tele- 
scope (UKIRT). The pixel scale is 0.281 arcsec pixel -1 which 
gives a field of view of ~72 arcsec. Our sample of quasars was 
observed during three observing runs in 09/1996, 09/1997 
and 05/1998. For the later two runs the image quality was 
exceptional with consistent FWHM of 0.45 arcsec observed. 

The K-band quasar images were taken using a quadrant 
jitter pattern. This cycled 2 or 4 times through a 4-point 
mosaic placing the quasar in each of the quadrants in turn. 
The actual position of the central value within each quad- 
rant was shifted slightly for each image to reduce the effect 
of bad pixels. Each image consists of ~ 100 sees of integra- 
tion time divided into exposures calculated to avoid satura- 
tion. The exposures varied between 5-10 sees for the quasars 
alone, down to 0.2 sees for the quasars with a bright star 
on the chip which we hoped to use as a psf star. Standard 
stars fr om the s ample of UKIRT faint standards (Casali & 
Hawarc en 1992 ) were observed for photometric calibration 
between observations of different quasars. All of the images 
were corrected for the non-linear response of IRCAM 3 using 
a formula supplied by the telescope support staff. 



3 OBTAINING THE CORRECT PSF 

Obtaining an accurate psf is vital to the analysis of the im- 
ages. With these ground-based observations the psf varies 
with seeing conditions and telescope pointing. An experi- 
mental psf was therefore determined for each of the quasars 
by observing a nearby bright star. This led to an unbiased, 
accurate psf without recourse to the quasar images. For three 
of the quasars, 0956-073, 1214+180 and 1216+069 there 
was a nearby star which could be placed on the frame with 
the quasar. This gave an accurate psf measurement with no 
loss of integration time on the quasar. If required, the posi- 
tion of the quasar for each observation was altered slightly to 
allow both the quasar and psf to be well within the bound- 
aries of the chip. For the remaining quasars the telescope was 
offset to a nearby bright star to use as the psf, before and 
after each quasar integration (which lasted a maximum of 
1600 sees). A number of psf measurements were therefore ob- 
tained for each night and each quasar. To ensure consistent 
adaptive optics correction, properties of the tip-tilt guiding 
were matched between quasar and psf measurements. To do 
this psf stars were selected to enable tip-tilt guiding from a 
star of a similar magnitude, distance from the object, and 
position angle to that used for the quasar image. Magni- 
tudes of the stars chosen to provide a psf measurement are 
given in Table |§| By examining fine resolution contour plots 
of the psf images, it was found that the psf was stable over 
the course of each night, but varied between nights at the 
telescope and for different telescope pointing. Because of 
this, the final stacking of psf images was performed with the 
same weighting between days as for the quasar images (see 
Section Q). 

As a test of the effectiveness of this procedure to provide 
the correct psf, the fit between measured psf and image for 
quasar 1543+489 has been compared to the fit between psfs 
measured for different quasars and the same image. Fig. [j] 
shows the radial profile of a 2 (image — psf) calculated in cir- 
cular annuli of width 0.5 arcsec. Here the psf has been scaled 
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radius/arcsec 

Figure 1. The radial profile of the variance (measured in annuli 
of width 0.5arcsec) of the difference between the image of quasar 
1 543+489 and the scaled psf obtained using the procedure in 
Section ^ (stars). For comparison, the profiles of the variance 
obtained using the f3 other psf measurements are also plotted 
(solid circles). The solid line is the best fit model to the variance, 
calculated as in Section 6.3 



so the total intensity is the same for both quasar and psf. As 
can be seen, the psf observed with the quasar image matches 
the quasar close to its centre better than any other psf. As 
the core of the psf is undersampled and the sampling be- 



tween psf and image has not been matched (see Section 
for further discussion of this), this result demonstrates the 
validity of the psf measuring technique. 



4 THE DATA REDUCTION 

The data reduction procedure was optimised to search for 
low surface brightness extended objects. The same proce- 
dure was used for both image and psf data so no extra dif- 
ferences between measured and actual psf were introduced. 
In calculating the flat-field for each mosaic, it was decided to 
ignore all pixels within the quadrant containing the quasar. 
This ensured that the flat-fielding technique was not biased 
to remove or curtail extended emission which could occur 
if a routine based on pixel values, such as a cr-clipping rou- 
tine, was used. Outside this quadrant, any areas occupied 
by bright stars were also removed from the calculation. The 
sky background level, assumed to be spatially constant was 
also calculated ignoring these areas. 

As the images were undersampled in the central regions 
we decided to use a sub-pixel shifting routine to centralise 
the images before they were median stacked to provide the 
final composite. Having replaced bad pixel values, the im- 
ages were shifted using a bicubic spline interpolation routine 
in order to equalise the intensity weighted centres, and were 
median stacked. Because the psf quality was found to vary 
from day-to-day, the final stacking of psf images was per- 
formed with the same weighting between days as for the 
quasar images. 

Finally, any nearby bright objects in the psf frame were 
replaced by the average in an annulus of width 0.5 arcsec at 
the distance of the object from the centre of the psf, around 
the centre. In the quasar frame any nearby objects were 



noted and blanked out of the error frame so they were not 
included when measuring \ 2 between image and model (see 
Section rat). 



5 DETERMINING SIMPLE MORPHOLOGY 

Because of the often disturbed morphology of quasar hosts it 
is not possible to immediately assume a form for the galaxy 
structure. For instance if the host is involved in a close 
merger, modelling it with a smooth profile will not provide 
the correct host luminosity. The extended wings of the psf 
from the intense nuclear component hide the host galaxy 
sufficiently that direct observation cannot easily reveal even 
violently disturbed morphologies. Simply subtracting a mul- 
tiple of the psf from the centre of the image will reveal some 
structure, but a deconvolution routine will reveal more struc- 
ture. The routine used was a modified Clean algorithm de- 
veloped for this problem which will be described elsewhere 
(Percival & Miller 200C). The results show that this routine 
was of sufficient quality to reveal the approximate symme- 
try of the host on a scale which includes most of the light 
important for modelling the galaxy. 

Examination of the deconvolved images revealed a clear 
distinction between disturbed and symmetric systems. Two 
of the quasars have morphologies which showed no sign of 
elliptical symmetry and instead show signs of recent merger. 
The deconvolved images of these quasars are shown in Fig. ^. 
From these images a value for the non-nuclear luminosity 
was obtained by summing the residual light excluding the 
central pixel. Unfortunately the non-nuclear structure re- 
vealed was not of sufficient quality to be extrapolated into 
the central region so the amount of nuclear light which orig- 
inates in the host galaxy is unknown. Magnitudes obtained 
from these deconvolved images should therefore be treated 
as approximate. Th e st ructure revealed for these quasars is 
discussed in Section 8.4. Deconvolving the remaining quasars 



revealed host galaxies with approximate elliptical symmetry. 



6 MODELLING THE QUASAR IMAGES 

Having determined that the extended structure around a 
quasar did not show signs of a disturbed morphology in- 
dicative of a close merger and revealed approximate ellip- 
tical symmetry, the luminosity and morphology of the host 
galaxy were estimated by fitting model images to the data. 
A x 2 minimisation technique described below was used to 
estimate the goodness of fit of the models. 



6.1 Producing a model galaxy 

In this Section we describe how the empirical galaxy surface 
brightness profile given by Equation |l| was used to estimate 
the contribution from the host to the counts in each pixel. 
This had to be done carefully because of the poor sampling 
of the images. The profile given by Equation ^ has proved to 



be an excellent fit to many d i fferent types of galaxy (Saon 



Capaccioli fc D'Onofrio 1993 ; Baggett, Baggett fe Anderson 



199c) and it is assumed that, if the hosts are not undergoing 



violent merger, this profile provides a good representation of 
the galaxy light. 
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Figure 2. The results of the deconvolution of quasars 1636+384 and 1700+518 revealing hosts with disturbed morphologies indicative 
of close merger events. The deconvolution output, obtained on the pixel scale, has been smoothed by convolving with a Gaussian with 
a = 0.5 pixels and the residual frame remaining after the algorithm has finished had been added back in to preserve luminosity. Contours 
are shown at 0.0125%, 0.025%, 0.05% and 50% of the peak intensity. For quasar 1700+581 a contour is also included at 0.00625% of the 
peak intensity. The size of the images is 7.3 X 7.3 arcsec (26 X 26 pixels). See Section %A for further discussion of the morphologies revealed. 



Before the method is described, it is useful to revise how 
an image is obtained from the light emitted by the quasar. 
Initially, the continuous distribution of light is altered by 
the atmosphere and the optics of the telescope in a way 
approximately equivalent to convolution with a continuous 
point spread function. The resulting continuous distribution 
is sampled by the detector which integrates the light over 
each pixel. This is equivalent to convolving the light with a 
square function of value 1 within a pixel and otherwise, 
and sampling the resulting distribution assuming uniform 
response across each pixel. The dithering and subsequent 
stacking of the images will provide another convolution, al- 
though by sub-pixel shifting the images prior to stacking, 
the effective smoothing width of this function is reduced to 
less than 1 pixel. The whole process can therefore be thought 
of as convolving the true psf, the quasar light and a narrow 
smoothing function (of width ~1 pixel) and sampling the 
resulting continuous image on the pixel scale. 

Because the psf measurements were obtained using ex- 
actly the same procedure as the quasar images, the measured 
psf is the result of a convolution of the true psf with the nar- 
row smoothing function. Convolution is commutative and 
associative so this smoothing function is accounted for in the 
measured psf and further smoothing of the model galaxy is 
not required. For this reason the unconvolved model galaxy 
should not be obtained by simply integrating the model pro- 
file over each pixel. Sampling the model galaxy profile onto a 
grid with spacing equivalent to the pixel scale and convolv- 
ing with the psf will not produce a correct model galaxy 
because of aliasing. 

In order to limit the aliased signal, the procedure 
adopted was to extrapolate the psf onto a grid which was 
finer than the pixel scale using a sine function (so no ex- 



tra high frequency components are introduced). The surface 
brightness of the model galaxy was then calculated at each 
point on a grid of the same size and was convolved with the 
psf on this grid. To provide the final model, this distribu- 
tion was subsampled onto the pixel scale. Progressively finer 
grids were used until the total counts in the sampled model 
galaxy converged, when the majority of the aliased signal is 
assumed to have been removed. The algorithm adopted used 
a fine grid with 4x the number of points at each successive 
step, and was stopped when the average of all the counts 
differs from that of the previous step by a factor less than 
0.01. 

Unless stated otherwise, all model host luminosities and 
magnitudes which relate to a 2D profile should be assumed 
to have been integrated to infinite radius. For the large ra- 
dius within which such models were fitted to the data, this 
makes only a small difference in the luminosity. The four pa- 
rameters of the host galaxy are the geometric radius of the 
elliptical annulus which contains half of the integrated light 
J?i/2, the total integrated host luminosity Lint, the projected 
angle on the sky a, the axial ratio a/b and the power law 
parameter (3. In this paper, the integrated host luminosity 
is quoted in counts /analogue data units (adu) detected in 
a 1 sec exposure. 



6.2 The nuclear component 

In principle, adding in the nuclear light is simple - the correct 
amount is added to the centre of the model galaxy to min- 
imise x 2 between model and observed images. However, it is 
important to account for all of the nuclear light. Differences 
between measured and true psf caused by undersampling, 
seeing variations or effects such as telescope shake must be 
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accounted for, even though the adopted observing strategy 
has limited some of thes e. In particular, when a nearby star 
(as used in Section 7.5) is deconvolved, the resulting light 
appears not only in the central pixel, but in the surround- 
ing pixels as well: if similar components from the nuclear 
light are not accounted for in the quasar images, the host 
luminosities and morphologies derived will be wrong. 

Because of the large peak in both the image and psf, 
trying to alter the sampling of the images by extrapolating 
onto a fine grid and resampling without inducing unwar- 
ranted frequency components causes 'ringing' in the images 
which is large enough to affect the results of the modelling. 
It would be possible to use a different extrapolation tech- 
nique, but this risks altering the image and measured psf 
in different ways. Instead, a more simple correction to this 
problem is adopted: rather than only adding a multiple of 
the psf to the central pixel, variable multiples of the psf 
are also added centred on the surrounding 8 pixels. In the 
perfect case where the measured psf is accurate, while such 
free parameters make convergence to the minimum \ 2 value 
slower, they do not affect the position of the minimum: 
the additional components make a negligible contribution 
to the model. However, suppose there is a discrepancy be- 
tween measured and true psf so that the deconvolved image 
of the nuclear light consists of a central spike surrounded 
by corrective components which decrease in magnitude with 
distance away from the spike. Allowing the value of the pix- 
els close to the core to be free parameters in our model will 
correct these discrepancies and any light observed originat- 
ing away from the core will be more likely to be from the 
host galaxy and not from escaping nuclear light. The oppo- 
site is also true, and these extra components will also correct 
psf measurement errors which cause light from the galaxy 
to be wrongly ascribed to the core (as for quasar 0956—073: 
see below). Undersampling problems do not affect the mod- 
elled galaxy to the same extent because the galaxy light is 
more uniformly distributed and discrepancies are smoothed. 
In particular, the total integrated light measured to be from 
the host will only be minimally affected: see below. 

Quasar 0956—073 was modelled using different numbers 
of these extra components, and the recovered host parame- 
ters are given in Table ^. As expected, the \ 2 value decreases 
with an increasing number of additive psf components show- 
ing that the fit between model and image is being improved. 
The recovered parameters for 1 or 9 additive components 
show moderate differences, but allowing more components 
makes no further significant change. Because the host lumi- 
nosity increases for this quasar with increasing numbers of 
added components, the sum of the extra psf contributions 
must be negative which suggests that, for this quasar, the 
psf has a slightly broader central profile than the quasar. 

For all of the quasars modelled, the total light within the 
eight extra psf components was not found to be systemati- 
cally positive or negative. If adding 8 extra 'psf components' 
around the core had always resulted in a total positive (or 
negative) component being subtracted from the quasar, this 
would have suggested that either these components were re- 
moving host light in addition to 'leaking' nuclear light, and 
that the host profile breaks down in a systematic way for 
these pixels, or that our observing strategy had produced a 
systematically incorrect psf. 

For all quasar host galaxy studies, there is no escap- 



psf 


Rl/2 




a/b 


a 


P 


x 2 


cmpts 


/kpc 


/adu 










1 


10.45 


312.7 


0.61 


1.05 


0.77 


2840.3 


9 


9.76 


350.1 


0.64 


1.04 


0.92 


2814.2 


25 


9.71 


353.3 


0.65 


1.03 


0.93 


2799.9 



Table 2. Recovered host galaxy parameters for quasar 0956—073 
modelled using different numbers of additive psf components. 



ing the fundamental problem that the galaxy profile has to 
extrapolated into the central region from some radius (to 
separate host and nuclear light). By adding in these extra 
components, all we're doing is extrapolating from different 
distances, and arguing that simply extrapolating only into 
the central pixel is not necessarily correct for these data. 
This is because the measured psf is incorrect for the (dis- 
crete) deconvolution problem we're trying to solve. 



6.3 The error frame 

Determining the fit between model and image requires an es- 
timate of the relative noise in each pixel, from both intrinsic 
noise in the image and differences between measured and 
true psf. Ideally these errors should be estimated without 
recourse to the images but, unfortunately, this is impracti- 
cal for these data. Faced with a similar problem, Taylor et al. 
( [1996D estimated the radial error profile by measuring the 



error in circular annuli of the image from which a multiple 
of the psf had been subtracted centred on the quasar with 
matched total luminosity. Using both the image and mea- 
sured psf in this way allows the error from psf differences to 
be included in the error frame. However, this model assumes 
that the host galaxies do not introduce any intrinsic varia- 
tions in the annuli within which the variance is calculated. 
Such variations could result from either differences between 
the radial profile of the host (convolved with the psf) and 
the psf profile, or significant deviation from circular host 
profiles. These effects will be small because the hosts only 
contribute a small percentage of the light and deviations 
from circular hosts are small. 

In order to reduce the number of parameters required to 
calculate the error profile and hopefully alleviate any dam- 
age caused by calcu lating the error frame from the data, 
Taylor et al. (1996) showed that a function of the form 
log(<Ti) = Aexp~ IJ - b<r/s) ~' -B, where A,B,S and 7 are four 
parameters, provides a good fit to the resulting profile. This 
profile models both the error in the central regions of the im- 
age and the Poisson background error outside the core. The 
four parameters are determined for each quasar by least- 
square fitting to the observed error profile. Such a fitting 
procedure also enables the error to be determined in the 
central regions where the gradient is too steep and there are 
too few points in each annulus to predict confidently the er- 
ror. In general this function fits the observed error profiles 
very well, and is used here (without including a contribution 
from the host) to estimate the errors in each pixel. 

Fig. ^ shows the observed and fitted error profiles for 
quasars 0956—073 and 1543+489 with and without includ- 
ing the best-fit model galaxy in the analysis. The best-fit 
host around quasar 1543+489 has an axial ratio of 0.89 in 
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the previous best-fit parameters until the total host lumi- 
nosity found for successive runs differed by less than O.f adu. 
The test ing performed for this algorithm is described in Sec- 
tion 







2 4 6 

radius/arcsec 



Figure 3. The calculated radial profile of the variance measured 
in annuli of width 0.3arcsec of the difference between the image 
of the quasar and the scaled psf (open trian gles ). The best- fit 
model for this profile calculated as in Section 3.3 is also plotted 



(dashed line). For comparison the radial profile, calculated in the 
same way, for the difference between the image and best-fit model 
galaxy with one nuclear component is plotted (solid triangles) 
and the model of this profile (solid line). This solid line shows 
the radial error profile used to produce the simulated data of 
Section The dotted line shows the Poisson noise level of the sky 
background. Top panel: quasar 0956—073, bottom panel: quasar 
1543+489. 



contrast to 0.64 for 0956—073. The deviation of the host 
around 0956—073 from circular symmetry explains why the 
error profile changes when including this host more than for 
quasar f543+489. 



6.4 x minimisation 

The algorithm used to find the global minimum in x 2 was a 
multi-dimensional direction set technique based on a method 
introduced by Powell in f964 ( Press et al. 199^ ). This algo- 
rithm requires an initial 'start point' from which it works its 
way downwards until it finds the minimum position. Briefly, 
the algorithm minimises \ 2 by sequentially adjusting each 
parameter (i.e. minimising along the axes of the parameter 
space), and then it minimises \ 2 on the vector along which 
the greatest change was made to x 2 m the previous steps. 
This procedure is repeated until the algorithm converges. 
Additionally, for all the quasars it was ascertained that the 
algorithm had found the correct minimum and not erro- 
neously finished early due to a numerical convergence prob- 
lem by repeatedly re-running the algorithm starting from 



7.2 



For all the results presented in Section ^, the algorithm 
was started from an initial position in parameter space corre- 
sponding to a broad, low luminosity galaxy. This was chosen 
so the algorithm avoided straying into a region of parameter 
space where all of the host light was in the core (i.e small 
R\/2)- This is a relatively flat region for \ 2 m parameter 
space and it can therefore take a long time for the algorithm 
to work its way out of this region. 

All pixels within a radius of 31 pixels (8.7arcsec) mea- 
sured from the centre of the quasar were included in the 
calculation of \ 2 ■ F° r au °f the best-fit host galaxies, the 
difference between model host luminosity within this area 
and the integrated luminosity was negligible, which implies 
that this area contained all of the important signal. 



7 TESTING THE MODELLING PROCEDURE 

7.1 Robustness to the error profile 

As a test of the robustness of the best-fit host luminosities 
to the determination of the error profile, we have modelled 
our image of quasar 0956—073 using different error profiles. 
Quasar 0956—073 was chosen for this test because the de- 
rived axial ratio of the host is the lowest of any qua sar (al- 
though the range of values is quite small: Section ^3). If the 
galaxy is important in the error frame calculatio n, th e error 

should 



profile calculated for this quasar as in Section 6.3 
be the most affected by the fact that we are ignoring the 
host (see Fig. ^|). Using an error profile calculated as in 
Section 6.3 but using the image only (i.e. not subtracting 
the psf), the integrated host brightness was found to drop 
from 350.1 adu to 316.6 adu, corresponding to a variation of 
^O.lmag. We have also tried re-calculating the error frame 
from the image minus the best-fit model image (galaxy and 
nuclear component convolved with the psf), again using the 
above formula to fit the error frame. Radial profiles of the 
two error frames are shown in Fig. ^. The best fit model 
parameters were used to calculate a new error frame, and 
we repeated this process until the best-fit host luminosity 
converged (subsequent iterations altered the integrated host 
luminosity by less than 0.1 adu). The final best-fit luminos- 
ity was found to be 351.5 adu, a negligible difference from 
the original minimum. 



7.2 Finding the minima 

Obvious tests to perform are that there is only one minimum 
for each quasar, and that the \ 2 function is well behaved 
around this point. Obviously, it is impossible to cover every 
position in parameter space to check that \ 2 is well-behaved 
and that there are no local minima. However, we have exam- 
ined the region of parameter space of interest using a variety 
of techniques and have found no potential problems. 

The minimisation algorithm is itself designed to cover a 
large region of parameter space; the algorithm sequentially 
searches for the minimum along a series of vectors (see Sec- 
tion 3.4 for details), and considers a large number of diverse 
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Figure 4. The variation of \ 2 (normalised to the minimum value) 
versus fixed total host luminosity. For each point on the curve, all 
parameters other than the luminosity have been altered to obtain 
the local minimum in x 2 ■ The 68.3% confidence interval, found 
by a separate binary search is also shown (dotted line). 



values along each vector. Rerunning the algorithm starting 
at the best-fit location previously found also tests any mini- 
mum along each axis in parameter space, as do es the calcu- 
lation of the error bars, described in Section 



7.4 



The shapes 

of the surfaces around each minimum are also revealed by 
this calculation. 

A test for local minima has been performed for quasar 
0956—073 over a larger region of parameter space: the min- 
imisation algorithm was started at a large number of di- 
verse initial host parameters, and no significant change in 
the best-fit parameters was discovered. Quasar 0956—073 
was chosen for this test because it has average signal-to- 
noise of any quasar. Fig. ^ shows a 'slice' through parameter 
space revealing how smoothly the constrained \ 2 minimum 
varies with fixed host luminosity for quasar 0956—073. To 
calculate each point of this curve, all of the parameters ex- 
cept the host luminosity were varied until the constrained 
X 2 minimum was reached. The remarkable smoothness of 
this curve demonstrates both that the global minimum is 
well pronounced and the function varies smoothly towards 
it, and that the minimisation routine is finding the correct 
minimum at each point: if it were not, a more rough surface 
would be expected, signifying that the optimum position 
had not been reached for each host luminosity. 



7.3 Using the \ statistic 

Use of the x 2 statistic is dependent on the error in each pixel 
being independent of the errors in the other pixels. This is 
expected if the errors in the images are dominated by Pois- 
son shot noise. Any large-scale differences between actual 
and model host could provide correlated errors, although 
these would hopefully have been discovered by the analysis 
of Section [|. It is possible that small-scale discrepancies re- 
main that extend across more than one pixel. However, the 
relatively large pixel scale works to our advantage by reduc- 
ing the likelihood of this. The central limit theorem then 



suggests that the error in each pixel should have approxi- 
mately Gaussian distribution. 

The minimum \ 2 values are highly dependent on the 
normalisation of the error frames, and cannot directly pro- 
vide tests of the model fits. The position of these minima 
are unaffected by the normalisation of the error frame as 
they are only dependent on relative variations between pix- 
els. Examining the reduced \ 2 values at the minima given 
in Table |i|, we see that the reduced \ 2 18 l ess than 1 for 
the majority of the quasars, and deduce that the procedure 
outlined in Section slightly over-estimates the error in 
each pixel. This is as expected due to the effect of t he h ost 
galaxy. The confidence intervals calculated in Section 7.4 will 



therefore be slightly too large, thus providing a moderately 
pessimistic error analysis. 

As any nearby companions were excluded when mea- 
suring \ 2 1 the number of pixels used, presented in Table ^, 
varies between quasars. For quasar 1214+180, a diffraction 
spike from a nearby star which ran close to the quasar was 
also excluded. Unfortunately the position of the pixels which 
were not modelled is more important than the number of 
such pixels and, for this quasar the position of the diffrac- 
tion spike was such that it covered a highly important region 
of pixels. Even though the area covered was small, the mod- 
elling suffered greatly. 

7.4 Calculating error bars on the parameters 

Provided that the galaxy model is a good representation 
of the true underlying host galaxy, the errors between the 
model and image are uncorrelated between separate pixels, 
and the procedure in Section |6~3| pro vides approximately the 
correct error frame (see Section |7.3|) , it is possible to calcu- 
late error bars on the true parameters using the x 2 statistic. 
The procedure to do this is to hold the chosen parameter 
fixed at a certain value, and minimise the remaining pa- 
rameters to find the local minimum in x' '■ The end points 
of the 68.3% confidence intervals on the best-fit parameter 
are given by the points for which A\ 2 = X 2 ~ Xmtn ~ 1> 
where Xmin ls the minimum value calculated al lowing all pa- 



rameters to vary (Bevington & Robinson 1992). A standard 
binary search has been used to find the required limits. As 
well as allowing error bars to be calculated, this procedure 
enables the parameter space to be examined and any prob- 
lems for each quasar to be spotted. 

In order to match the light from the host galaxies, the 
behaviou r of the integrated host luminosity, j3 a nd -R1/2 are 
coupled (Abraham, Crawford & McHardy 1992). The deter- 



mination of the error bars is therefore complicated by the 
question 'what limits, if any should be placed on the pa- 
rameters being adjusted to find the constrained minima?'. 
In finding the global minima, all of the parameters are ef- 
fectively allowed to vary over all space: although bounds 
are placed on the parameters, they are not reached (except 
when modelling the star, see Section ^^). However, at fixed 
integrated host luminosity, these limits are often reached be- 
cause the profiles required to optimally match the light do 
not necessarily have to be those of galaxies. The philosophy 
adopted is that all the parameters should be allowed to vary 
except f3, upon which limits of 0.25 < f3 < 6.0 should be set 
to provide some adherence to standard galaxy profiles. 

For quasar 0956—073, we have examined the required 
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cut through parameter space for the integrated host lumi- 
nosity, calculated by minimising all other parameters to ob- 
tain each point. The distribution of local \ minima are 
shown in Fig. ^: the curve displays simple structure, mono- 
tonically decreasing to the global minimum from both direc- 
tions so we are justified in using the simple A\ 2 = 1 cut-off 
for the error bars. The resulting 68.3% confidence interval 
for the luminosity is also shown. 

The value of x 2 depends on the error frame used, and 
it is expected that the error bars do so as well. The effect 
of altering the error frame for quasar 0956—073 has been 
tested by using the error frame calculated from the image 
only as in Section 3.3. Using this error frame, the 68.3% 
confidence interval on the host magnitude changed slightly 
from -25.11 < Mk < -24.87 to -25.10 < M K < -24.66. 



7.5 Fitting to a normal star 

On the same frame as quasar 0043+039 we observed a star 
of similar signal-to-noise as the quasar. As a test of the fit- 
ting procedure we decided to see if we could fit a 'galaxy' to 
the star. Starting from an initial position in parameter space 
corresp onding to a broad, low luminooity galaxy aa adopted 



in all of the modelling, the resulting best-fit parameters arc 

given in Table ^. As can be seen, the fitting procedure rolled 
down the hill towards a host galaxy of very low luminosity. 
At such low total luminosity, the remaining four galaxy pa- 
rameters are poorly determined: altering these parameters 
results in a very small change in \ 2 ■ Consequently it is no 
surprise to find that the best-fit = 6 value is one of the 
limits set in the modelling procedure. 



8 RESULTS OF THE ANALYSIS 
8.1 Luminosities 

For three of the quasars, analysis of how x' 2 varies within 
the parameter space revealed that the best-fit host luminos- 
ity was not well constrained. A host galaxy was determined 
as being present in that a lower limit was determined in all 
cases. However, the maximum light which could have come 
from the host was not clear because the shape of the host 
was not sufficiently resolved. The morphology of the best- 
fit galaxy at large Lint could alter to place the majority of 
the host light in the central region. This effect could have 
been avoided by placing limits on R1/2 or, for instance, using 



the near-infrared Fundamental Plane (Pahre, Djorgovski & 



de Carvalho 1998), although these upper limits would have 



been highly dependent on the criteria set. The host luminos- 
ity is ultimately limited by the total light in the image, and 
it is expected that the host luminosities for these quasars 
do have upper bounds at high values of Lint, but these high 
values would not be of any use in determining the actual 
host light. 

For the remaining nine quasars, the minima were suf- 
ficiently constrained to provide 68.3% confidence intervals. 
Comparison of different confidence intervals provided infor- 
mation on the depth of the valleys within which each mini- 



compare with the LLband host galaxy studies undertaken by 
McLeod & Rieke (1994a;b), we convert their total (nuclear + 
host) and host luminosities to the Jf-band by applying a sin- 
gle conversion factor to the apparent magnitudes. This then 
sets the relative normalisation of the Tf-band and LLband 
samples; conversion to absolute magnitudes is subsequently 
undertaken in exactly the same way for all of the infra-red 
samples. 

In a study of the energy distribution of the PC quasars 
(from which McLe od & Rieke chose their samples), Neuge- 
bauer et al. ( |l987|) foun d (H - K) = 0.90 for the sample of 
Mcleo d fc Rie ke ( 1994a| ) and (H - K) = 0.98 for McLeod & 
Rieke ( f!994b| ). In the upper panel of Fig. [|, we adopt these 
values to convert the total luminosities of the McLeod & 
Rieke quasars into the Lf-band. 

The light from the galaxy component is assumed to be 
dominated by an evolved stellar population, the colour of 
which reddens with increasing redshift. For nearby galax ies, 
H-K ~ 0.25, which was used by McLeod & Rieke (09951) to 



convert galaxy absolute magnitudes. However, the apparent 
H — K is dependent on redshift and, at the redshifts of the 
quasars imaged by McLeod & Rieke (1994a; b), H — K ~ 0.6 



is expected for an evolved stellar population ( Lilly fc Longair 



1984). This was adopted to convert the McLeod & Rieke 
galaxy luminosities into the K-haxiA. 

We have also checked the calibration of the McLeod & 
Rieke samples and our sample of modelled quasars (with 6 
overla pping objects) against the data of Neugebauer et al. 
(1987). The average total quasar luminosity for the subsam- 



ples are in good agreement, although individual values vary 
by up to 0.7 mag, presumably due to intrinsic quasar vari- 
ability. 

On e quasar (1354+213) was i maged by McLeod & Rieke 



(1994b), Neu gebau er et al. (1987) and in our study. Neuge- 
bauer et al. (1987) derived H — K — 1.0 for this object, 
which is higher than H — K = 0.3 derived by combining 
the McLeod & Rieke LLband and our if-band observation. 
However, the McLeod & Rieke and our observations were 
undertaken at different epochs, and the luminosity is not 
expected to remain constant. 



The study of Taylor et al. (1996) was performed in the 



Lf-band, and the apparent A'-band magnitudes of host and 
nuclear components were taken directly from this work. The 
data from the different infra-red samples were then con- 
verted to absolute m agnit udes, applying the /^-correction 
of Glazebrook et al. ( |l095| ) for the host galaxy and assum- 
ing the nuclear component follows a standard power law 
spectrum f(u) — v~ ' 5 . 



Using the error bars calculated in Section 7.4 to weight 



the data, the average integrated host galaxy magnitude for 
our quasars was found to be (Mk) = —25.15+0.04. For com- 
parison, when converted for c osmo logy exactly as our data, 
the sample of Tay lor et al. ( [l996|) gives (M K ) = -25.68, 
McLeo d fc R ieke ( |l994a| ) (Mk) = -25.42 and McLeod & 
Rieke ( |l994b| ) (M K ) = -25.68. 

Recent determinations of the K-band luminosity of an 
L* galaxy (Gardner et al. 1997) have resulted in M K = 
— 24.6, compared to previo us determinations of My _ —24.3 



mum was found and the quality of each determination. 



(iGlazebrook et al. 1991 ) and M K = -25 .1 (|Mobasher 



Host and nuclear luminosities for our quasars are com- 
pared with the results of other studies in Fig. [| In order to 



Sharpies fc Ellis 1993| )"~The Gardner et al. ( |1997| ) valuels 
plotted in the top panel of Fig. |^. This shows that the av- 
erage luminosity of our hosts is ~ 1.6 times that of an L* 
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quasar -R1/2 ^int axial ratio a (3 nuc/host ^fhost Ktot K ps { 

/kpc /adu /radians ratio 



0043+039 


7.41 


355.0 


0.96 


1.44 


0.60 


13.9 


16.04 


13.54 


11.90 


0137-010 


5.09 


407.5 


0.66 


2.55 


2.05 


13.7 


15.96 


13.46 


10.83 


0244-012 


3.28 


165.6 


0.86 


2.66 


0.61 


16.6 


16.87 


14.18 


11.15 


0316-346 


8.56 


864.8 


0.76 


0.07 


1.25 


12.2 


15.15 


12.73 


9.62 


0956-073 


9.76 


350.1 


0.64 


1.04 


0.92 


21.4 


16.02 


13.07 


9.20 


1214+180 


3.24 


259.7 


0.74 


1.25 


1.12 


15.2 


16.34 


13.75 


8.52 


1216+069 


19.75 


379.3 


0.74 


2.42 


2.73 


21.4 


15.93 


12.99 


11.10 


1354+213 


11.55 


481.4 


0.69 


1.38 


0.73 


6.79 


15.67 


13.82 


9.21 


1543+489 


7.75 


320.1 


0.89 


3.00 


0.67 


25.9 


16.12 


12.98 


9.92 


1636+384 












24.1 


17.65 


14.11 


10.07 


1700+518 












14.8 


15.86 


11.82 


10.23 


2112+059 


5.59 


169.1 


0.93 


0.49 


2.13 


67.0 


16.85 


12.72 


11.08 


2233+134 


8.01 


175.0 


0.80 


1.54 


0.64 


28.5 


16.88 


13.64 


10.17 


2245+004 


5.16 


572.9 


0.81 


1.93 


3.16 


1.25 


15.52 


14.87 


10.02 


star 


16.6 


0.000031 


0.16 


2.26 


6.0 











Table 3. Best-fit host galaxy parameters as determined by the 2D modelling described in Section Q Also included for comparison are 
the best-fit parameters for a nearby star found on the frame of quasar 0043+039 which had similar signal to noise as the quasar. For 
quasars 1636+384 and 1700+518 deconvolution of the images revealed a highly disturbed morphology which extended close to the core 
of the quasar and model fitting was not attempted. Consequently host galaxy magnitudes .presented for these quasars are the relatively 
inaccurate measurements calculated from the deconvolved images as described in Section tj. Nuclear-to-host ratios are calculated in the 
rest frame of the quasar from the derived absolute magnitudes in order that these values are consistent with Fig. ^| (see Section 8.1 for 
details). The apparent magnitudes of the quasar, the host component, and the star used to give a psf measurement are also presented. 



quasar \ 2 number of reduced Li nt /adu (host) Mx(tot) 







pixels modellec 


x 2 


min 


best-fit 


max 


min 


best-fit 


max 






0043+039 


2639 


79 


2885 





92 


287.7 


355.0 


454.4 


-25 


56 


-25 


29 


-25 


07 


-28 


22 


0137-010 


2661 


25 


2987 





89 


333.7 


407.5 








-25 


08 


-24 


87 


-28 


01 


0244-012 


2985 


64 


2974 


1 


00 


63.9 


165.6 


999.4 


-26 


85 


-24 


90 


-23 


86 


-28 


01 


0316-346 


2674 


08 


2987 





90 


766.3 


864.8 


1009.7 


-25 


61 


-25 


41 


-25 


31 


-28 


24 


0956-073 


2814 


17 


2987 





94 


315.8 


350.1 


394.4 


-25 


11 


-24 


98 


-24 


87 


-28 


36 


1214+180 


2499 


93 


2509 


1 


00 


112.7 


259.7 








-24 


94 


-24 


03 


-27 


96 


1216+069 


2437 


36 


2625 





93 


290.2 


379.3 


604.0 


-25 


62 


-25 


11 


-24 


82 


-28 


49 


1354+213 


2281 


11 


2987 





76 


462.4 


481.4 


502.1 


-25 


20 


-25 


16 


-25 


11 


-27 


38 


1543+489 


2567 


80 


2874 





89 


263.6 


320.1 


402.3 


-25 


56 


-25 


31 


-25 


10 


-28 


89 


2112+059 


2876 


00 


2945 





98 


37.7 


169.1 








-24 


91 


-23 


29 


-29 


50 


2233+134 


2679 


73 


2941 





91 


108.7 


175.0 


269.2 


-24 


58 


-24 


11 


-23 


59 


-27 


78 


2245+004 


2584 


98 


2987 





87 


438.1 


572.9 


941.6 


-26 


24 


-25 


70 


-25 


41 


-26 


48 


star 


3055 


42 


2969 


1 


03 




0.000031 





















Table 4. Table showing the end points of the 68.3% confidence intervals calculated as in Section [7,4| for the best-fit integrated host 
luminosities. To further aid comparison between model fits for different quasars, the \ 2 values for the best fit model are presented with 
the number of pixels used in this calculation. Note that the position of the un-modelled pixels is more important than the number of 
such pixels. 



galaxy. Note that for all three values, the derived average 
luminosity is 1 — 2 time s that of an L* galaxy, and the con- 
clusions of Section |ll. l| are not affected by this choice. 

We compare our sample to recent HST i?-band results 
in the lower panel of Fig. |H| assuming an apparent R — K = 
2.5 for the total light from our quasars based on the average 
value for the 6 quasars whic h ove rlap our sample and the 
sample of Neugebauer et al. (1987). The R — K colour of an 



evolved stellar population, assumed to dominate the host 
galaxies, is dependent on the redshift of the source and, for 
the re dshifts of our samp le (z ~ 0.35), is expected to be 
~ 3.5 ( Dunlop et al. 1989[ ). All the data (including our data 
after conversion to apparent i?-band magnitudes) presented 
in the bottom panel of Fig. ^ were adjusted for cosmology 
assuming that the nuclear component has a spectrum of the 



form f(y) oc v 



and the galaxy component has f(y) oc 



8.2 Morphologies 

Morphologies are parametrised by the best-fit value of (5: 
13 — 1 values correspond to disk-like, and f3 = 4 to spheroidal 
profiles. The technique described in Section 7.4 has been 



used to reveal how well the (3 parameter is constrained by 
the modelling. The result of this analysis is presented in Ta- 
ble |. As can be seen, the fi parameter is well constrained for 
fewer quasars than the luminosity and \ 2 error bars reveal a 
highly skewed distribution for the expected true value given 
the best-fit f3 value. In order to correctly determine the dif- 
ferential probability between disk and spheroidal profiles, we 



© 0000 RAS, MNRAS 000, 000-000 



Quasar host galaxies 11 



-24 



M K (nuclear) 
-26 -28 



-30 



CNI 

I 

CO 
CNI 
I 



* CM 



Csl 
I 



CM 
I 

to 

CM 
I 



a: cm 
S CM 



CM 
I 



A 


+ 


+ 


+ 






x^ 






++•- 


a"4 x 


.+ 4 



^ 



□ , 



4U J -4 



-22 



-24 
M R (nuclear) 



-26 



Figure 5. The top panel shows nuclear vs. int egrated host absolute i\-band magnitudes for our sample of quasars (solid circles) with 
error bars calculated as described in Section The errors in the measured nuclear component are derived from these and consequently 
the err ors ar e strongly correlated. Plotted for co mpariso n are the calculated host magnit udes fo r the radio-quiet quasars imaged by Taylor 
et al. (1996) (open triangles), McLeod & Rieke ( |l994a| ) (crosses) and McLeod & Rieke (|1994q) (plus symbols). Details of the conversion 



The luminosity of an L* galaxy, Mi 



-24.6 



of the McLeod Sz Rie ke data from the /f-band to the A'-band can be found in Section S.l 
(Gardner et al. 1997) is also plotted (dashed line), as is the locus of points with a rest-frame K"-band nuclear-to-host ra tio of 8 (dotted 
line). In order to compare with recent HST results, in the bottom panel we convert our data into the R-band (see Secti on for details'). 
Symbols for our data are as for the top panel. The ij-band best-fit luminosities from disk or exponential galaxies (Hooper, Impey & 
Foltz 19p^) are also plotted (diamonds separated by dotted lines). No at tempts were made to distinguish the best host morphology in 
this work. The derived ij-band host luminosities of McLure et al. (1999 ) are also sho wn (radio-quiet sample: open squares, radio- loud 
sample: solid triangles). The luminosity of an L* galaxy, = —21.8 (Lin et al. 1996) is plotted for comparison (dashed line). 



need to know the relative dispersion of (5 for each morpho- 
logical type. However, examining the best-fit parameters, 
the error bars on (3, and the shape of x 2 surface, on which 
we have information from the binary search to find the error 
bars, we can infer the best fit morphology for some of the 
quasars. The suggestion from this is that luminous radio- 
quiet quasars can exist in hosts dominated by either disk- 
like or spheroidal components. A histogram of these data is 
plotted in Fig. ^| where the distribution is compared to that 
recovered from simulated data with exact (3 = 1 or f3 = 4 
profiles. 



8.3 Axial ratios and angles 

Analysis of the parameter space reveals that the axial ratio 
and projected angle of each host are better constrained than 
the other parameters. Fig. ^| shows histograms of these pa- 
rameters for all quasars modelled. The distribution of axial 



ratios is small with (a/6) = 0.79 ± 0.03. This is in agree- 



ment with those found by McLure et al. ( 1999), but higher 
than found by Hooper, Impey & Foltz (1997). The projected 



angles are uniformly distributed as expected. 



8.4 Highlighted results for selected quasars 

8.4.I Quasar 0043+039 

The broad-absorption-line (BAL) quasar PG0043+039 has 
been subject to 2 previous studies to determine host galaxy 
propertie s. It w as observed in the i band by Veron-Cetty & 
Woltjer dl990[ ) who determined M % = -23.9 if the host is 
a disk like Jj3 = 1) galaxy, or Mi = —24.7 for a spheroidal 
(P — 4) galaxy. This quasar was also ob served using the 
wide- field camera on HST by Boyce et al. (1998), who used 
a cross correlation technique to determine that the host was 
slightly better fit by a disk galaxy with My ~ —21.6. We 
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quasar 


min /3 







max /3 


morphology? 


0043+039 


0.42 





60 


0.86 


disk 


0137-010 


1.52 


2 


05 


>6.0 


spheroid 


0244-012 


<0.25 





61 


1.62 


disk 


0316-346 


1.02 


1 


25 


3.62 


? 


0956-073 


0.71 





92 


1.21 


disk 


1214+180 


<0.25 


1 


12 


>6.0 


? 


1216+069 


1.44 


2 


73 


>6.0 


spheroid 


1354+213 


0.65 





73 


0.82 


disk 


1543+489 


0.49 





67 


0.90 


disk 


2112+059 


<0.25 


2 


13 


>6.0 


? 


2233+134 


<0.25 





64 


1.21 


disk 


2245+004 


2.17 


3 


16 


5.55 


spheroid 



Table 5. Table showing 68.3% confidence intervals for the best- 
fit host f) parameters for the 12 quasars modelled using the 2D 
\ 2 minimising techniq ue (S ection H). Error bars were calculated 
as described in Section 7.4 with (3 constrained to lie in the range 



0.25 < /3 < 6.0. The morphology of the dominant contribution to 
the host galaxy is also presented, based on the best fit f3 parameter 
and associated confidence interval. 



0.2 0.4 0.6 0. 
Axial ratio 




12 3 

Projected angle/radians 

Figure 6. Histograms showing the distribution of axial ratios 
and projected angles of the host galaxies. These parameters 
were well constrained for all of the 12 quasars modelled, and 
are therefore plotted from all of the minima found. Axial ra- 
tios are tightly constrained with with a/b > 0.64 for all hosts, 
and (a/b) = 0.79 ± 0.03. The distribution of projected angles is 
approximately uniformly distributed. 



also find that the dominant morphology is disk-like and cal- 
culate M k — -25.29. The old burst model of Bruzual & 
Chariot (1993) predicts V — K = 3.3 which is consistent 
with the derived V - K = 3.7. 



8.4.2 Quasar 0316- 346 

This quasar was previously observ ed using the wide field 
camera on HST (Bahcall et al. 1997) and the host was found 
to reveal evidence of a merger, in particular tidal ta ils ex - 
tending ~ 20kpc west of the quasar. Bahcall et al. (1997) 
also provide a 2D fit to the host properties and find that the 
best-fit host is disk galaxy with My — —22.3. We also cal- 
culate a best-fit disk galaxy and find Mk = —25.44, giving 
V — K = 3.1, again consi stent with the old burst model of 
Bruzual & Chariot (|l993[) which predicts V - K = 3.3. 



8.4.3 Quasar 1214+180 

There have been no previous attempts to determine the mor- 
phology of the host galaxy around this quasar possibly due 
to the nearby star which was utilised in this work to ob- 
tain an accurate psf. Unfortunately in our images, a diffrac- 
tion spike from this star passed close to the quasar reducing 
the area that could be used to calculate \ 2 - Although the 
modelling converged to give basic galaxy parameters, further 
analysis of the parameter space revealed that this minimum 
was not well constrained. 



8.4.4 Quasar 1216+069 

Our analysis of this quasar benefited because the images 
were obtained using the tip-tilt system and there is a nearby 
bright star which was placed on the same frame as the 
quasar and used to obtain a psf measurement. Pre viously, 



'nebulosity' has been observed arou nd this quasar (Hutch 



ings, Crampton & Campbell 1984), and a more detailed 
HST study fou nd a best-fit sphero idal (/3 = 4) galaxy with 



M v = -22.3 ( [Boyce et al. 1998| ). We also find that the 
most likely host is a large spheroidal galaxy and obtain 
M K = -25.1, giving V - K = 2.8. 



8.4.5 Quasar 1354+213 



Using a psf subtraction technique, McLeod & Rieke (19941) 
found a residual host galaxy with Mk = —25.6 when 
converted to our cosm ology using the /("-correction from 
Glazebrook et al. (1995) a nd t he apparent colour correction 
H-K = 0.6 (see Section |tJ). Our best fit host luminosity 
was Mk = —25.2. Analysis shows that the luminosity and 
the 13 parameter are both tightly constrained by the mod- 
elling and the best-fit (3 = 0.73 suggests that the host is dom- 
inated by a disk component. The rest-frame nuclear-to-host 
ratio for this quasar is only 6.8 (the apparent nuclear-to-host 
ratio is 4.6), which explains why the derived parameters have 
small error bars. 



8.4.6 Quasar 1636+384 

We are not aware of any previous attempts to determine the 
luminosity and morphology for the host galaxy of quasar 
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1636+384. Preliminary deconvolution of the light revealed 
that the excess, non-central light displayed a morphology 
greatly disturbed from elliptical symmetry (as shown in 
Fig. 0). The structure includes an excess of light to the NW 
of the core which is interpreted as a merging component as 
well as light around the central core which probably origi- 
nates from the host. From this image it was unclear how to 
distinguish between the host and the interacting companion, 
so the luminosity of the host was estimated by summing 
pixel values excluding the central pixel. This provided an 
approximate if-band absolute magnitude of Mk = —23.5. 



8.4.7 Quasar 1700+518 

Quasar 1700+518 is a bright BAL quasar of low redshift 
(z = 0.29). Such low redshift BAL objects are rare and hard 
to discover since the broad absorption lines are in the UV 
and consequently quasar 1700+518 has received much inter- 
est: specific studies of this quasar have been undertaken in 
many different wavebands ( Hutchings, Ncff fc Gowcr 1992 



Stockton, Canalizo & Close 19981; Fines et al. 19991). Be 



cause of the low redshift and the brightness of the quasar, 
1700+518 has also been included in many sampl es of quasars 
imaged to obtain details of their host galaxies (Neugebauei 



et al. 1985 



McLeod 



Rieke 1994b| ), although these have 
only provided upper limits for the host magnitude. More 
recent imaging studies have shown that the morphology of 
the underlying structure consists of a disturbed host pre- 
dominantly to the SW of the core (Stockton, Canalizo & 



Close 1998) and a close interac ting companion to the NE 
flHutchings, Ncff fc Gowcr 1992| ) which is most likely a ring 



galaxy ( |HinT!?r*(T*a7"T99*S ) . Deconvolution of the light from 
this quasar, as shown in Fig. |i| confirms this picture of the 
structure. With the disturbed morphology it is difficult to 
know how to split the light in the central pixels into nuclear 
and host components. As for quasar 1636+384 the host lu- 
minosity was estimated by summing the counts in the pix- 
els surrounding the central one (ignoring those from the NE 
companion) . There will be errors caused by leakage of light 
from the nuclear component and from the contribution of 
the host to the central pixel. An approximate Tf-band ab- 
solute magnitude of —24.9 was obtained for the host galaxy 
and —24.4 for the NE companion galaxy. 



8.4.8 Quasar 2233+134 



Both Smith et al. (1986) and Veron-Cetty & Woltjer (199C) 



included this quasar in their samples, but both failed to re- 
solve the host galaxy beyond o btainin g upper limits for the 
luminosity. Hutchings & Neff (1992) did resolve the host 



galaxy and found the host to be best-fit by a (3 — 4 model, 
although they did not resolve further information about the 
galaxy. However, we find that the most probable host has an 
disk profile and calculate Mk = —24.1, the lowest luminos- 
ity host modelled. If we constrain the host to have an ellip- 
tical profile, the best fit luminosity becomes Mk = —25.9, 
although the half light radius is very small for this model 
{R1/2 = 1.5 kpc) which places it a long way from the K- 
band Funda mental Plane of Pahre, Djorgovski & de Car- 
valho (1998). If the host parameters are constrained to lie 
on this plane, then rerunning the modelling gives a best fit 



host with Mk = —24.6. Neither of these changes would be 
sufficient to alter our conclusions. 



9 SIMULATED DATA I - SINGLE 
COMPONENT GALAXIES 

Trying to recover known host parameters from Monte-Carlo 
simulations of the actual data enables the distribution of re- 
covered parameters given the true values to be determined. 
Note that the error bars calculated using the x 2 statistic 
are instead determined from the distribution of possible true 
values given the data. These two distributions are not nec- 
essarily equal. We need to determine the distribution of re- 
covered values in order to answer questions such as 'Are our 
results biased towards low (3 values?'. 

In view of the distribution of recovered (3 values, it was 
decided to simulate data to match the images of quasars 
0956—073 and 1543+489. These quasars span the distribu- 
tion of signal-to-noise of all the images, and 2D model fitting 
revealed evidence for disk dominated hosts in both cases. 
Verification of this result is interesting as recent work has 
suggested that the hosts of luminous quasars shou ld be dom- 
inated by the spheroidal component (see Section 11.2). 



9.1 Creating the mock data 

Simulated galaxi es w ere created using the procedure out- 
lined in Section 6.1 and a single 8 function added to the 
centre of each to create a 'perfect unconvolved model'. The 
height of the 8 function was chosen to match the total sig- 
nal of the original images. These models were then convolved 
with the psf measured to match the quasar. 

Gaussian noise was added with a radially depend ent 



variance as given by the error profile calculated in Section 6.3 
including the best-fit host galaxy in the calculation. The er- 
ror profiles used for quasars 0956—073 and 1543+489 are 
given by the solid lines in Fig. ||. Differences between mea- 
sured and true psf were included in this analysis, and are 
therefore included in the noise levels added to the simulated 
data. This noise model assumes that the errors in different 



pixels are independent (see Section 7.3). 

We have simulated 100 images with exact disk hosts, 
and 100 images with exact spheroidal hosts for each of the 
two quasars chosen. The true integrated host luminosity was 
set at 300 adu for simulated data of both quasars. This con- 
servative value is below the best-fit value obtained from the 
data for both quasars, providing a stringent test of the mod- 
elling. This is particularly true for a f3 = 4 host: constraining 
P — 4 when modelling the observed image would have re- 
sulted in a best-fit Lint ~S> 300 adu. The simulated images 
were analysed using exactly the same 2D modelling proce- 
dure described above for the observed data. The range of 
recovered parameters is analysed below. 

9.2 Results from the simulated data: luminosities 

Recovered luminosities, presented in Fig. ^ reveal a skewed 
distribution, particularly for hosts with exact spheroidal 
profiles where the recovered luminosity is biased towards 
a low value. This is consistent with the morphology being 
skewed towards a low j3 value (see next Section): if (3 is 
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quasar j3 Li llt /adu 

true min average max true min average max 



0956 


-073 


1.0 


0.74 


1.14 


2.28 


300.0 


262.1 


323.7 


510.0 


0956 


-073 


1.0 


1.93 


3.77 


7.42 


300.0 


213.6 


288.3 


463.1 


1543- 


f489 


1.0 


0.79 


1.08 


1.69 


300.0 


260.1 


324.7 


512.2 


1543- 


f489 


4.0 


2.13 


4.04 


7.25 


300.0 


209.6 


302.2 


469.8 



Table 6. Table showing the mean and 68.3% confidence intervals for the recovered fi and integrated host luminosity from the simulated 
data. 100 simulations were performed for each morphology for each quasar. 
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Figure 7. The distribution of recovered luminosities from Monte-Carlo simulations of 100 /3 = 1 images and 100 f} = 4 images. The j/-axis 
gives the number of recovered values within each luminosity bin. Noise has been added to match ob serv ations of quasar (a) 0956—073 and 
(b) 1543+489, including a contribution from the error in the measured psf as described in Section The luminosity of each simulated 
galaxy, marked by the dotted line, was set at 300 adu. 



decreased, the luminosity also has to decrease to keep the 
counts in the outer pixels (those most important for fitting 
the host) the same. The counts in the centre of the galaxy 
are less important because of the additional nuclear compo- 
nent which is adjusted to match the data. 

The mean and variance in the recovered luminosities 
are presented in Table ^. Although the error bars reveal the 
extent of the skewed distribution, the mean is within 10% 
of the true value for each quasar and morphology. 

9.3 Results from the simulated data: 
morphologies 

The skewed distribution observed in the error bars on the 
true host f3 value is mirrored by the distribution of f3 values 
recovered using the standard modelling procedure described 
in Section^. Fig. ^ shows the relative distribution of (3 values 
retrieved from the simulated images. Limits of 0.25 < f3 < 8 
were placed on fitted /3 values. For quasar 0956—073, 16 of 
the simulated images created with exact spheroidal hosts, 
had recovered ft > 8. For quasar 1543+489, this number was 



14: these values are not included in Fig. |8j The distribution 
was used to calculate the mean and standard deviation given 
in Table ^, assuming all fits with f3 > 8 actually had j3 = 8. 

If the host were a spheroidal galaxy with (5 = 4, the 
probability of recovering a best-fit value of /3 < 1 is ~ 0.03 
for 0956-073 and ~ 0.01 for 1543+489: the best-fit values 
from the images were f3 = 0.92 and /3 = 0.67 respectively. 
The evidence for the existence of hosts dominated by a disk 
component therefore appears to be strong. In Fig. ^| the dis- 
tribution of retrieved /3 values for the 12 quasars modelled 
is also shown. This distribution is inconsistent with the hy- 
pothesis that all the hosts are dominated by spheroidal com- 
ponents on the scales probed by these measurements. The 
histogram is divided to show the probable distribution of 
morphologies given the options (3 = 1 or (3 — 4. As can be 
seen, the modelling suggests that approximately half of the 
hosts are dominated by disk components. 



© 0000 RAS, MNRAS 000, 000-000 



Quasar host galaxies 15 



(a) 0956-073 



(b) 1543+489 



(c) all quasars 



o 



o 

CNI 



o 



o 

CNI 



model 0=1 
(disk) 



o 



o 

CM 



model 0=4 
(spheroid) 



o 



o 

CM 



~\ — 1 — r 



model 0=1 
(disk) 



EbL_ 



J , L 




model 0=4 
(spheroid) 



r>r>unrk 



2 4 6 
best fit 

[_| probable disk 
probable spheroid 



4 




4 





Figure 8. Evidence that the hosts of luminous radio-quiet quasars are not exclusively dominated by spheroidal components. The 
distribution of f3 values recovered from 2-D modelling of 100 simulated images created with hosts using exact disk or spheroidal profiles 
is presented. Noise has been added to these images to match observati ons of (a) quasar 0956—073 and (b) quasar 1543+489, including a 
contribution from the error in the measured psf as described in Section |s),l| . The dotted line shows the best-fit /3 value recovered from the 
images of these quasars, (c) For comparison the distribution of best-fit /3 values obtained from all of our i\-band images is also plotted. 
The probable morphology of the host was determined from the x 2 error bars derived for the true value given the data. 



10 SIMULATED DATA II - TWO 
COMPONENT GALAXIES 

In order to constrain the potential importance of spheroidal 
cores in the galaxies found to be dominated by disk-like 
profiles, we have analysed synthetic quasars created with 
two host galaxy components. Using the Fundamental-Plane 
(FP) relation betw een R1/2 and Lint found in the K-band by 
Pahre et al. (1995), we have added extra spheroidal (/3 = 4) 
components to the recovered best-fit host galaxy of quasar 
1543+489. Note that this best fit host had f3 = 0.67. We have 
tried the same analysis using ft — 1 and found no change in 
the effects produced by the spheroidal core. After adding in 
the nuclear component and noise as described in Section ^| 
we have recovered the best-fit host galaxy parameters us- 
ing our single component modelling. Spheroidal components 
were added with a variety of different luminosities, and five 
different realisations of the additional noise component were 
added to each. The resulting average recovered Lint & /3 are 
given in Table |^. 

Because R1/2 (spheroidal) and Li„t (spheroidal) follow 
a FP relation, the importance of this component is en- 
hanced for large Li n t (spheroidal) and diminished for small 
Lint (spheroidal). The recovered total luminosity for small 
spheroidal components is therefore very similar to that of 
the disk alone. For large spheroidal components, the mod- 
elling places an excess of host light in the core in order to 
simultaneously fit the outer disk-like profile and the inner 
profile with a single, large /3 value. This explains the be- 
haviour of the difference between the actual and recovered 
Lint values. Recovered f3 monotonically increases with the 
increasing luminosity of the spheroidal core, suggesting that 





Mnt 


/adu 




Rl/2 





spheroidal 


total 


recovered 


diff 


/kpc 




0.0 


320.1 


291.2 


-28.9 


8.26 


0.61 


40.0 


360.1 


306.5 


-53.6 


8.19 


0.63 


80.0 


400.1 


337.9 


-62.2 


8.02 


0.69 


120.0 


440.1 


382.7 


-57.4 


7.71 


0.81 


160.0 


480.1 


436.3 


-43.8 


7.35 


0.96 


320.0 


640.1 


785.8 


145.7 


5.45 


1.93 


480.0 


800.1 


1305.2 


505.1 


3.85 


3.37 



Table 7. Average recovered host parameters from single com- 
ponent fits to synthetic quasars with 2-component host galaxies. 
Uncorrclated Gaussian noise has been added to these models to 
match that of quasar 1543+489, and the average recovered val- 
ues are given for 5 different realisations of this noise. The same 
noise was added to corresponding mock images created with dif- 
ferent spheroidal luminosities, so there will be a systematic error 
because 5 realisations are not sufficient to fully sample the recov- 
ered parameters with the given noise level. The result of analysing 
a host with no spheroidal component shows that the results sys- 
tematically underestimate /3 and Li n t and overestimate R\/2 by 
small amounts. Note that the relative dependence of the recovered 
parameters on the spheroidal component will not be affected. 



the spheroidal core cannot be completely 'hidden' without 
affecting the best fit galaxy. This adds to the evidence that 
the low f3 values recovered for some of the quasars implies 
that they do not contain strong spheroidal components. Note 
that the recovered host luminosities are approximately cor- 
rect for recovered values of j3 consistent with a host domi- 
nated by a disk-like profile. 
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For the quasars which have best-fit hosts dominated 
by spheroidal components, a disk-like profile at larger radii 
could have erroneously increased the recovered total host 
luminosity. However, in order to simultaneously fit these re- 
gions, small values of R1/2 were required. For the quasars 
with hosts found to be dominated by spheroidal components, 
the relatively large values of R1/2 recovered suggest that 
such a disk-like component is not present. 



11 DISCUSSION 
11.1 Luminosities 

The integrated host luminosities derived from our Tf-band 
images exhibit a low dispersion around a mean similar to 
that calculated in studies of less luminous quasa rs. This is 
in accord with the work of McLure et al. (1999) who also 
found no evidence for an increasing trend, although they 
had fewer data points at high nuclear luminosity. 

Previous HST studies have found evidence that host 



luminosity increa ses with nuclear luminosity (Hooper, Im 



pey & Foltz 1997), although the trend observed in this work 



could be due to incorrect nuclear component removal: escap- 
ing nuclear light which increases in luminosity with the core 
could be added to the host light. It has recently been stated 
that the psf derived by pack ages s uch as tinytim, as used 
by Hooper, Impey & Foltz (1997) deviate from empirical 
WFPC 2 psfs at lar ge radii (> 2 arcsec ), due to scattering 



within the camera ( McLure et al. 199!: ), and this could be 
the reason for an excess of nuclear light at larger radii which 
could be mistaken as host light. This excess light could also 
be the reason th e low axial ratios observed in the Hooper, 
Impey & Foltz ( 1997 ) w ork a re not in accord with those 
derived in McLure et al. or in this A'-band study. 

The triangular shape of the McLeod & Rieke points in 
Fig. |E] found for low redshift (0 < z < 0.3) Seyferts and 
quasars of lower luminosity than those in our sample, has 
been shown to be in accord with a lower limit to the host lu- 



minosity whic h increases with nuclear luminosity (McLeod 



& Rieke 1995). This cut-off in the host luminosity is equiva- 
lent to there being an upper limit to allowed nuclear-to-host 
ratios. The triangular shape is not followed by the results 
of the work presented in this paper which lie to the right 
of the McLeod & Rieke points. The relative positions of the 
two data sets in this Figure are set by the empirical H — K 
corrections applied to the apparent //-band data (see Sec- 
tion i.l for details). Quasar 1354+213 was included in bot h 
our sample and the sample of McLeod & Rieke (1994b), 



and the results of both studies independently suggest a rest- 
frame nuclear-to-host ratio of 7 — 9. This places 1354+213 
at the right of the triangular shape of the McLeod & Rieke 
points in Fig. ||, but it has a nuclear-to-host ratio lower than 
most of the quasars in our sample, and is therefore to the left 
of most of our points. We co nclude that the limit suggested 
by McLeod & Rieke ( 1995 ) must break down for quasars 
with the highest nuclear luminosities. 

This is in con trast to recent work by McLeod, Rieke & 
Storrie-Lombardi ( 1999 ) who claim that the lower bound on 
host luminosity extends to the highest luminosity quasars, 
partly based on the discovery of one luminous quasar, 
1821+643 which appears to be in a host at ~ 25 L* . What 



should we expect? The hosts of the quasars known to date 
already extend to about 2L* . Should the hosts of quasars 
which are ten times more luminous be found in galaxies at 
20L*? Our analysis suggests not. 

This result is highly important for recent quasar m odels . 
In particular, the model of Kauffmann & Haehnelt (199£) 



predicts that the upper limit to the nuclear-to-host ratio 
should extend to quasars such as those imaged in this work. 
However, this is clearly not the case. A possible fix to their 
model would be to invoke the scatter of the Magorrian et al. 
relations to explain high luminosity quasars (& high 



([l99§) 



mass black holes) within low luminosity structures, and in- 
voke a steeply-declining host mass function to explain the 
lack of really massive hosts. Further work on this model 
would then be required, particularly with regard to the re- 
vised slope of the high luminosity tail of the quasar lumi- 
nosity function. Alternatively, factors other than black-hole 
mass, such as nuclear obscuration, accretion processes, etc. 
could be the cause of differing nuclear luminosities within 
reasonably similar galaxies (with similar black hole masses) . 

11.2 Morphologies 

Recent HST results suggest that luminous nuclear emission 
predomi nantly arises from h osts with large spheroidal com- 
ponents ( McLure et al. 1999| ) . The two least luminous radio- 
quiet quasars imaged by McLure et al. (199£) have disk- 



like structure at radii > 3 arcsec, while the more luminous 
quasars are completely dominated by spheroidal profiles. 
Could we be seeing a relationship between host morphol- 
ogy and nuclear luminosity? This is particularly interesting 
when compared to the black hole mass-spheroid mass and 
spheroid mass-luminosity rel ations determined for nearby 
galaxies by Magorrian et al. ( |199S ) : a large black hole, po- 
tentially capable of powering luminous AGN appears more 
likely to be present in galaxies with large s phero idal com- 
ponents. Both the results of McLure et al. (199E) and the 
relations of Magorrian et al. (1998) suggest that quasars 
with strong nuclear emission should predominantly exist in 
hosts with large spheroidal components which dominate any 
disk-like structure. 

By careful analysis we have provided evidence that 
a large fraction of the host galaxies found in this work 
are dominated by disk-like profiles. However, the most im- 
portant light for this modelling comes from radii greater 
than those of the HST study, where the disk component, if 
present, is expected to be strong. The if-band images de- 
scribed here are not of sufficient quality for us to resolve 
the inner region and produce a 2-component fit to the host 
galaxy. This is in contrast to results from HST where the 
increased resolution enables the inner region to be resolved, 
and the spheroidal core of the galaxy becomes more impor- 
tant for modelling with a single j3 parameter. By analysing 
synthetic data, we have been able to show that for hosts 
where we find a dominant disk-like component, any addi- 
tional spheroidal component will not result in a large change 
in the recovered total host luminosities. We have also pro- 
vided suggestive evidence that the spheroidal cores of these 
quasars are of relatively low luminosity. Further analysis of 
both the regions and profiles probed by different studies, 
and higher resolution data on the cores of the quasars anal- 
ysed in this work would be very interesting, and could help 
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to explain the different morphological results of recent host 
galaxy studies. 



12 CONCLUSIONS 

We have presented the results from a deep K-band imag- 



ing study designed to reveal the host galaxies of quasars 



with higher luminosities than targeted by previous studies. 
Extending host-galaxy studies to these quasars was made 
possible by the stability provided by the tip-tilt adaptive 
optics system at UKIRT, which enabled accurate psf mea- 
surements to be obtained for the deep quasar images. We 
have been able to resolve host galaxies for all of our sample. 

The principle conclusion of this study is that the lumi- 
nous quasars in this sample have host galaxies with similar 
luminosities to quasars of lower total luminosity. Derived 
nuclear-to-host ratios are therefore larger than those ob- 
tained by previous work, and place these quasars beyond 
the upper limit suggested by studies of quasars with lower 
total luminosities. Host morphologies are less certain, but 
there is weak evidence that the hosts of these quasars can 
be dominated by either disk-like or spheroidal profiles on 
the scales probed by these images. 
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